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The evaporation of a drop from a liquid subphase is investigated. The two liquids 
are immiscible, and the contact angles between them are given by the Neumann 
construction. The evaporation of the drop gives rise to flows in both liquids, which 
are coupled by the continuity of velocity and shear-stress conditions. We derive 
self-similar solutions to the velocity fields in both liquids close to the three-phase 
contact line, where the drop geometry can be approximated by a wedge. We focus 
on the case where Marangoni stresses are negligible, for which the flow field consists 
of three contributions: flow driven by the evaporative flux from the drop surface, 
flow induced by the receding motion of the contact line, and an eigenmode flow that 
satisfies the homogeneous boundary conditions. The eigenmode flow is asymptotically 
subdominant for all contact angles. The moving contact-line flow dominates when 
the angle between the liquid drop and the horizontal surface of the liquid subphase is 
smaller than 90°, while the evaporative-flux driven flow dominates for larger angles. 
A parametric study is performed to show how the velocity fields in the two liquids 
depend on the contact angles between the liquids and their viscosity ratio. 
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I. INTRODUCTION 



The evaporation of a liquid drop from a flat solid substrate has been studied exten- 
sively, for example in the context of evaporation-driven particle deposition 1-1 and evapora- 
tive cooling . However, apparently little is known about the evaporation of liquid drops 
from a liquid subphase. On a liquid subphase the contact line of the drop is not pinned, 
the subphase is deformable, and the evaporation can generate a flow in both liquid phases. 
Examples of such a system are gasoline or oil drops on water, or water drops on mer- 
cury. Similar interface deformations occur for drops evaporating from soft, deformable solid 
substrates 8-11 , for which the flow geometry is reminiscent of drops floating on another liquid 
subphase. 

Evaporation of sufficiently small sessile drops under atmospheric conditions is usually 
governed by the diffusive transport of vapour in the surrounding air 1 ' 12-17 . The evaporative 
flux from the drop surface drives a radially outward flow inside the drop, which is responsible 
for the so-called coffee-stain effect 1 ' 2 . In drops where the contact line is free to move, 
there is an additional contribution to the flow field that comes from the receding motion 
of the contact line 14,18 . Furthermore, the non- uniform evaporation of sessile drops leads to 
temperature gradients on the liquid-air interface. The subsequent Marangoni flows can give 
rise to additional circulation 5 ' 19 ' 20 . 

It has remained a challenge to characterize the flow near the contact line, due to the small 
length scales involved and the evaporative singularity. Recently, we have studied the nature 
of the flow near the pinned contact line of a liquid drop evaporating from a flat solid surface 
by approximating the drop profile by a wedge geometry 21 . We have shown that the Stokes 
flow in this wedge-shaped region can be described by similarity solutions and consists of 
three contributions: flow driven by the evaporative flux from the drop surface, flow induced 
by the downward motion of the liquid-air interface, and an eigenmode flow that satisfies the 
homogeneous boundary conditions and can give rise to Moffatt corner eddies 22 . 

Here, we investigate the evaporation of a liquid drop from a liquid subphase, as shown 
in Fig. 1(a). We consider drops with a size smaller than the capillary length. In that case, 
the interface of the liquid subphase remains horizontal and the drop attains the shape of a 
lens. The angles 9, a, and /3 between the liquids are dictated by the three surface tensions 
in the system via the Neumann construction, such that both the horizontal and vertical 
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FIG. 1. (a) A drop of liquid 1 evaporating from a non-volatile liquid subphase (liquid 2). The 
angles 9, a, and f3 between the two liquids can be found from the Neumann construction at the 
three-phase contact line. The dashed square marks the area close to the contact line, where the 
drop shape can be approximated by a wedge, (b) An overview of the wedge geometry containing 
both liquid 1 and liquid 2. The contact line is located at the origin of the polar coordinate system 
(p, (J)). The evaporative flux J (indicated by the arrows) from the surface of liquid 1 drives the 
flow inside liquid 1, which is coupled to the flow in liquid 2 via the continuity of velocity and shear 
stress conditions. 

components of the capillary forces add to zero 23 . We focus again on the nature of the flow 
near the contact line of the drop, where the drop shape can be approximated by a wedge 
geometry, see Fig. 1(b). The evaporation of liquid 1 will cause a flow in both liquids 1 and 
2. We describe these flows by deriving self-similar solutions to the Stokes equations in a 
wedge geometry, in which the flows in both phases are coupled dynamically. Again, the flow 
in the wedge consists of three contributions: the evaporative-flux driven flow, the flow due 
to the receding motion of the contact line, and the eigenmode flow for the coupled phases. 
The eigenmode flow in a similar double-wedge region was already investigated by Anderson 
& Davis 21 . They found conditions under which Moffatt eddies can be present in viscous 
flow of two fluids near a corner of two rigid planes. The present work is a variation of this 
problem, where there are two free surfaces and one of the liquids evaporates. A parametric 
study is performed to investigate how much flow is generated in the liquid subsphase by the 
evaporation of the liquid drop as a function of the angles 8, a, and /3, and the viscosity ratio 
of the two liquids. We first treat the three flow contributions separately, and then identify 
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the regime in which each contribution is asymptotically dominant. 



II. PROBLEM FORMULATION 

The geometry of the drop close to the contact line can be approximated by a two- 
dimensional wedge, see Fig. 1(b). The three contact angles between the liquids in the 
wedge add to 9 + a + f3 = 2tt. To ensure both horizontal and vertical balances between the 
surface tensions at the three-phase contact line, these angles are restricted to 9, a, (3 < tt . 
We adopt a polar coordinate system (p, 0) with the origin located at the contact line and 
define the velocities in terms of streamfunctions *i(p, 0) and *2(p ; 0) as 

lMi <9*i 19*2 9*2 

The flow is governed by the Stokes equations, or equivalently, in terms of the streamfunctions, 
by the biharmonic equations 

V 4 *! = 0, V 4 * 2 = 0. (2) 
The general self-similar solutions to the biharmonic equations in the two liquids are 25 

*i(p, 4>) = p m [ai cos m0 + a2 sin m0 + 03 cos(m — 2)0 + 04 sin(m — 2)0] , (3a) 
* 2 (p?0) = p n [c\ cosn0 + c 2 sinn0 + c 3 cos(n — 2)0 + c 4 sin(n — 2)0] , (3b) 

where m,n 7^ 0,1,2; in that case, the form (3) is degenerate and additional terms arise. The 
exponents m and n will be selected by the boundary conditions. 

The solution to the flow field described by (2) consists of three contributions 21 . One 
contribution comes from the evaporative mass flux from the free surface of liquid 1, which 
drives a flow inside the drop. A second contribution arises from the volume change of the 
drop due to the evaporative mass loss. In case the contact line of the drop is pinned, there 
will be a corresponding contribution to the flow field due to the downward motion of the 
liquid-air interface 21 . On a liquid subphase however, the contact line of the drop is not 
pinned. Instead, there is a receding motion of the contact line, which results in a decrease 
in the drop radius and gives a contribution to the flow field. The third contribution to 
the flow field is the classical "eigenmode" solution to the homogeneous problem for which 
all boundary conditions are zero, as described by Dean & Montagnon 26 , Moffatt 22 , and 
Anderson & Davis 21 . From now on, we will refer to these contributions as the flux solution, 



the moving contact-line solution, and the eigenmode solution. Since (2) is a linear system, we 
can consider each of these three contributions separately. The full solution can be obtained 
by superposition. 



A. Flux solution 

For the evaporative-flux driven flow, we impose as boundary conditions that liquids 1 
and 2 are immiscible and the geometry is fixed. At the interface between liquids 1 and 2, 
located at = (or (j) = 2n) this gives 



Uh{p,0) = 0, 
m 02 (p,2tt) = 0. 



(4a) 
(4b) 



The liquid- air interfaces, located at <fi = 9 (liquid 1) and <p = 9 + a (liquid 2), are shear-stress 
free 



T P <t> x {p,Q) =0, 



(5a) 
(5b) 



The coupling between liquids 1 and 2, at the boundary = 0, is made by the continuity of 
velocity and shear stress conditions at the liquid-liquid interface, 



or 7/i 
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T p(l>2 (p, 2tt), 
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V2 
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P2 



p dcf) 



(6a) 



(6b) 



=2tt 



with 7/i and 7/ 2 the dynamic viscosities of liquids 1 and 2, respectively. As a consequence, 
the viscosity ratio 7/ = 7/1/7/2 will appear as a parameter in the solution. The conditions 
above are the homogeneous boundary conditions for the coupled system of \l/i and ^2- The 
inhomogeneous condition comes from the evaporative flux that drives the flow. We impose 
that liquid 1 evaporates, while liquid 2 is assumed non-volatile 



u^(p,0) 
Ufoip^a + 6) 



1 



J(P), 



(7a) 
(7b) 



with J the evaporative flux and pe the liquid density. The evaporative flux from the drop sur- 
face is known from the solution to the diffusion equation describing the vapor concentration 
field in a corner geometry 1 



The prefactor A can be obtained from the asymptotics of the full spherical-cap solution 
(Popov 13 ), and A(a) = n/2a. Here U = DAc/Rpt is the velocity scale, which is of order 
pm/s for water drops under atmospheric conditions', with R the radius of the drop as 
indicated in Fig. 1(a), D the diffusion constant for vapor in air, and Ac = c s — the 
vapor concentration difference (in kg/m 3 ) between the drop surface and the surroundings. 
Hence, for a > ir/2 the evaporative flux diverges as the contact line is approached, while for 
a < it/2 the evaporation is suppressed near the contact line. 

The eight boundary conditions for the biharmonic equations (2) are thus given by (4-7). 
The inhomogeneous evaporative flux condition (7a) is driving the flow in liquid 1, which in 
turn generates a flow in liquid 2 via the coupling conditions (6). From condition (7a) we 
find that exponent in (3a) is m = X(a). The exponent in (3b) is selected by the coupling 
condition (6a), which implies n = m = A(a). The coefficients a,\ to C4 in (3a) and (3b) can 
be found from the boundary conditions (4)-(7). The system of equations that has to be 
solved to find these coefficients is 



J{P) 



A(a)U p^- 1 . 



(8) 



Cx = b, 



(9) 



with 



x = 



(di, 0,2, CI3, CL4, Ci, C2, C3, C4) , 



(10) 



a vector containing the coefficients a± to C4. The right-hand side vector b contains only 
zeros except for element 6 7 = AU, which comes from the evaporative driving term. The 
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coefficient matrix C is given by 



C = 
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(11) 



A solution to (9) can be found when det(C) 7^ 0. Critical points in the solution to (9) appear 
for those values of a, 9, 77 where det(C) = 0. In these cases the eight boundary conditions 
are no longer linearly independent and the eigenmode solution comes into play 21 , as will be 
discussed in Section II C. 

The 8x8 matrix given by (11) describes the coupled flow in the two liquids. We anticipate 
that a decoupling into two 4x4 systems arises for rj = and for 77—7-00. In these limiting 
cases, the flows in the two liquids are completely independent. The boundary conditions 
experienced by each of the two liquids in these decoupled cases are sketched in Fig. 2. For 
?7 = 0, liquid 1 effectively experiences a no-slip condition at the interface with liquid 2, as 
if it were contacting a solid substrate. This feature is due to the extremely high viscosity 
of the liquid subphase. For 77—7-00, the stresses in liquid 2 are so small that effectively a 
no-stress condition applies and the interface is a free surface. For liquid 2 it is exactly the 
other way around: no-stress at the interface with liquid 1 when 77 = and no-slip when 
77—7-00. 
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FIG. 2. Schematic of the boundary conditions experienced by liquids 1 and 2 in the two limiting 
cases, (a) r\ = 0: Liquid 1 experiences a no-stress condition at the free surface and a no-slip 
condition at the interface with liquid 2, whereas liquid 2 experiences two no-stress conditions, (b) 
r] — > oo: Liquid 1 experiences two no-stress conditions whereas liquid 2 obeys a no-stress condition 
at the free surface and a no-slip condition at the interface with liquid 1. 

B. Moving contact-line solution 

As the drop evaporates its volume decreases and the contact line recedes. The resulting 
translation of the interface of liquid 1 gives rise to kinematic boundary conditions at the 
interface between liquids 1 and 2 (located at = 0) and the free surface of liquid 1 (located 
at = 9) 



with R(t) the radius of the drop as indicated in Fig. 1(a), and dR/dt the contact-line 
speed. We solve the biharmonic equations describing the flows in both liquids subject to the 
boundary conditions (13) and (5-6), (7b). Since (13a) and (13b) are independent of p, we 
obtain for the exponents in (3) m = n = 1. For m,n = 1, the solutions (3a) and (3b) are 
degenerate, and the solutions to the biharmonic equations are of the form 22 




(13b) 



(13a) 



(p, 0) = p [a\ cos + a 2 sin + a 3 </> cos (ft + a 4 0sineV] , 
^2 (p, 4>) — P [ c i cos + c 2 sin + c 3 cos + c 4 sin 0] . 
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(14a) 
(14b) 



By solving for the coefficients in (14a) and (14b) using boundary conditions (13) and (5-6), 
(7b), we obtain 

dR 

^(p, <P) = tr 2 (p, <P) = —p sin(a + 6-(j>). (15) 

This result is simply a uniform flow away from the contact line. 

To find an expression for the contact-line speed dR/dt we can again use the solution to 
the evaporative outer problem 13 

dV j,DAc 

— = -txR /(«), (16) 

at p£ 

with V the drop volume. An expression for / can be obtained from the integrated evaporative 
flux from a spherical cap and is given by lb 

. sin(7r - a) A [°° 1 + cosh2(7r - a)r . . . . 

f(a) = - —, ^ + 4 / — ^ '— tanh ar dr. 17 

v ' l + cos(vr-a) J smh27rr v ' K ' 

In contrast to the model of Popov 13 , the drop geometry on a liquid subphase consists of two 
spherical caps instead of one. In this case, the drop volume is given by 



V(a,8) = 7iR 3 



-cos 3 (a) + 3cosfa) + 2 cos 3 (a + 9) - 3 cos(a + 9) - 2 



3 sin a 3 sin (a + 



ttR h(a, 9). 

(18) 



Assuming that contact angles a and 9 remain constant as the contact line recedes, we obtain 
an equation for the rate of change of the drop radius 

^ = -luJ^L (19) 

dt 3 h(a,ey 1 ; 

where U = DAc/Rpi is again the characteristic velocity, which is known for a given drop 
geometry. From (19) we obtain a prediction for the evolution of the drop radius in time that 
could be checked experimentally 



with tf the total evaporation time of the drop. 



C. Eigenmode solution 



The eigenmode solution is the nontrivial solution to the homogeneous problem 22 . The 
form of the solution is still given by (3a) and (3b), but the exponents m and n are selected 



by the criterion det(C) = 0, with C given by (11). In case the determinant of the coefficient 
matrix C is equal to zero, the system of equations is degenerate and one of the equations 
becomes redundant. In practise, one can find the coefficients of the eigenmode flow by 
setting one coefficient equal to an arbitrary value (here taken unity) and removing one line 
from the coefficient matrix (11) to remove the redundancy. The last column of the coefficient 
matrix then serves as the right-hand side vector. This procedure implies that the strength 
of the eigenmode flow cannot be determined from the "inner" Stokes flow problem close to 
the contact line, but is set by the outer flow in the drop 21,27 . 

For the standard case of a drop evaporating from a flat solid substrate, the condition for 
the (Moffatt) eigenmode solution was found to be 21 

M(X E , 9) = sin 2(X E - 1)6 - (X E - 1) sin 26 = 0, (21) 

where X E is the exponent of the eigenmode solution. To ensure regularity of the velocity 
field at the origin, only the roots that have Re(X E ) > 1 are considered. The case A^ = 2 is a 
trivial root and does not correspond to an eigenmode. The present case of a drop evaporating 
from a liquid subphase is more general, since the substrate is a liquid that is not flat. The 
limit of r]=0 effectively corresponds to a liquid evaporating from a solid substrate, and in 
this case det(C) reduces to 

det(C)=F(a,\ E )M(\ E ,e), (22) 

where 

F(a,0) = -8Asin{A[27r- (a + 0)]} sin {A [2?r - (a + 6)} + 2(a + 6)} . (23) 

When a = n — 6, F = — 8Asin 2 (ttX) and we exactly recover the case of a liquid evaporating 
from a flat solid substrate 21 given by (21). For arbitrary 77 and a the equation det(C) = 
gives the solution to the generalized corner-flow problem of a liquid on a liquid subphase, 
where, instead of a no-slip condition, a coupling condition between the two liquids applies. 
Hence, the exponent of the eigenmode solution is a function of the three problem parameters: 
m = n = X E (a, 6, 77). This problem is a variation of the one studied by Anderson & Davis 24 , 
where a two-dimensional viscous flow of two fluids in the corner between two solid walls was 
studied. Here, instead of two solid walls there are two free surfaces. 
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III. RESULTS 



A. Flux solution 

The flows in the evaporating drop and the liquid subphase depend on the angles between 
the phases a, 9, and /3, and the viscosity ratio 77. To illustrate the effect of these parameters 
on the velocity fields in both liquids, we show some representative cases. To start we 
extend our previous work where we studied the flow inside a liquid evaporating from a solid 
substrate 21 . The current problem is more general, but will tend to the flat solid case when 
^<1, such that liquid 1 experiences a no-slip condition at the interface with liquid 2, and 
a = 7r — 9, such that the interface of liquid 2 is flat. In Fig. 3(a) we show that in this limit 
we indeed recover the same flow patterns as described in 21 , with the flow in liquid 1 being 
much stronger than the flow in the more viscous liquid 2. Similarly to the results obtained 
in our previous work, we observe that for larger contact angles 9 (and hence smaller a) the 
flow shows a sequence of reversal structures. When 9 = 50° (a = 130°), the flow is separated 
into two regions, such that near the interface with liquid 2, the flow in liquid 1 is directed 
away from the contact line, towards the center of the drop. The dashed line is the separating 
flow line that ends in a stagnation point at the contact line. The smaller the angle a, the 
more separatrices appear. Very little flow is generated in the liquid subphase, which is very 
viscous in the limit rj <C 1. The influence of the liquid subphase becomes apparent when 
the viscosity ratio rj is increased, see Fig. 3(b). When the viscosity of liquid 2 becomes 
comparable to or smaller than the viscosity of liquid 1, the evaporation of liquid 1 generates 
a flow in liquid 2 via the coupling conditions (6a) and (6b). The change in viscosity ratio 
also affects the flow in liquid 1, which now experiences a different boundary condition at 
the liquid-liquid interface. This results in a disappearance of the separatrices in liquid 1 for 
9 = 50°, 128°, and 130°. 

Note that the cases shown in Fig. 3 are somewhat artificial, since a liquid subphase would 
not remain flat: the Neumann condition at the contact line will require that a + 9 > 180°, to 
ensure a vertical force balance. To explore a physically more realistic case, we now consider 
the situation where all surface tensions are equal, for which 9 = a = (3 = 120°. In Fig. 4 we 
plot the resulting streamline patterns for three viscosity ratios rj. It can be observed that 
the velocity field in liquid 1 nearly obeys a no-slip condition for r/ = 1/10, which means that 
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FIG. 3. (a) Streamline plot of the flux solutions *$>i (dark gray, blue online) and ^2 (light gray, 
red online) for a viscosity ratio r/ = 1/10 (liquid 2 much more viscous than liquid 1). Contact 
angle 9 is varied from 50° to 160°, whereas a is adjusted such that the interface of liquid 2 is fiat: 
a = 180° — 9. The dashed lines represent the separatrices between the different flow regimes, (b) 
Same as in (a) but with a viscosity ratio rj = 10, which means that liquid 1 is much more viscous 
than liquid 2. 



there is almost no evaporation-driven flow in liquid 2. When the viscosity ratio increases, 
the evaporation of liquid 1 generates more flow in liquid 2, which can be seen from the larger 
arrow size in Fig. 4 for r] = 10 compared to 77 = 1/10. For r] = 10, the flow in liquid 1 follows 
an almost free-slip condition. The strength of the flow generated in the liquid subphase as 
a function of the viscosity ratio will be analyzed further below. 

In Fig. 4 we do not see the reversing flow structures we observed in Fig. 3. The reason 
for this feature is that the sequences of wedge-like flow structures shown in Fig. 3 are 
triggered by large exponents A. We can understand this behavior by inspecting the self- 
similar expressions for the streamfunctions given by (3): since m = n = A, a high frequency 
in the 0-direction arises when A>1. Hence, when a -C 90°, which means A > 1, multiple 
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FIG. 4. Streamline plot of the flux solutions (dark gray, blue online) and ^ (light gray, red 
online) for = a = (3 = 120°, which means all surface tensions are equal. The viscosity ratio is 
varied from r] = 1/10 (liquid 1 experiences almost no slip at the interface with liquid 2) to rj = 10 
(almost free slip). 




FIG. 5. Streamline plot of the flux solutions (dark gray, blue online) and \&2 (light gray, red 
online) for a = 60° and a = 20°, where 6 is varied from 165° to 175°, for r/ = 10. For small a, a 
sequence of wedge- like flow structures is observed in both phases. 

wedges appear. To illustrate this behavior we show the flow structures for large and small 
a in Fig. 5. For a = 60°, only one separatrix is present in the flow in liquid 1. When a is 
decreased to 20° however, we obtain a sequence of four wedges, which implies that the flow 
is frequently changing directions. 

We have seen that the evaporation of liquid 1 generates a flow in liquid 2 through bound- 
ary conditions (6a) and (6b). The strength of the flow generated in liquid 2 depends on 
the viscosity ratio, as shown in Fig. 3 and 4. A measure for the amount of flow generated 
in liquid 2 by the evaporation of liquid 1 is the ratio u P2 (<f) = a + 8)/u ( / )1 ((f) = 9). This 
ratio relates the amount of flow generated at the free surface of liquid 2 to the driving term, 
which is the flow towards the free surface of liquid 1, see the inset in Fig. 6(a). Fig. 6 shows 
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FIG. 6. (a) Plot of the flow generated at the free surface of liquid 2 compared to the driving by 
the evaporative flux, u^/ui, as a function of the viscosity ratio. The inset shows the definition of 
ui = u<f>i(<ft = Q), the velocity towards the free surface of liquid 1, and U2 = u P2 (cj) = a + 9), the 
velocity along the free surface of liquid 2. Here, a = 120°, whereas 9 is varied from 100° (black, 
solid line) to 120° (gray, dashed line) to 140° (light gray, dotted line), (b) Same velocity ratio as 
in (a) but now as a function of 9. Parameter rj is fixed to 1, whereas a is varied from 100° (black, 
solid line), to 120° (gray, dashed line) to 140° (light gray, dotted line). 

how this velocity ratio depends on the parameters r], 9, and a. In Fig. 6(a) we observe 
that when rj = 0, the flow in liquid 1 obeys a no-slip condition at the interface between the 
two liquids, which means no flow is generated in liquid 2. In the opposite limit r] — > oo a 
free-slip condition applies, with a maximum amount of flow generated in liquid 2. The flow 
strength in liquid 2 depends on the viscosity ratio between the two liquids, and also on the 
contact angles 9 and a. In Fig. 6(b) it is shown that not only the magnitude but also the 
flow direction depends on 9 and a: M2/M1 can change sign. This sign change is due to the 
appearance of separatrices in the velocity field, as observed in Fig. 5. 

B. Moving contact-line solution 

The moving contact-line solution is just a uniform flow away from the contact line, as 
given by (15). The corresponding flow pattern is shown in Fig. 7. No gradients in the 
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FIG. 7. Streamline plot of the moving contact-line solutions \&i (dark gray, blue online) and ^ 
(light gray, red online) for various a, 6, and a viscosity ratio rj = 1. 

velocity fields in the two liquids are observed, which means there is no viscous stress. By 
contrast, the hinge contribution to the velocity field in a pinned evaporating drop on a solid 
substrate leads to a significant viscous stress, due to the pinning of the contact line and the 
no-slip condition at the liquid-solid interface- 1 . As the exponents of the streamfunctions 
for the moving contact-line solution are constant, the flow patterns shown in Fig. 7 are 
the same for all angles a and 9 and viscosity ratios 77. While the exponent of the moving 
contact-line solution is equal to unity, the exponent A for the flux solution is larger than 
unity when a < tt/2. Therefore, we anticipate that the moving contact-line solution will be 
the dominant motion near the contact line when a < tt/2. 

C. Eigenmode solution 

The streamline patterns of the eigenmode solution are shown in Fig. 8. In the plots, 
we took rj — 1, such that the eigenmodes in liquids 1 and 2 are coupled via the boundary 
conditions (6a) and (6b). We observe that the velocity field in liquid 1 is directed inwards 
along the interface with liquid 2, and outwards along the free surface. The same flow 
structure applies to liquid 2, however, for a = 20° a reversal in the flow direction in liquid 2 
is observed. In that case, the liquid-liquid interface acts as a separatrix. For the parameters 
values used in the plots in Fig. 8 the exponents of the eigenmode solutions are real-valued, 
which means no viscous eddies are present 22 . 

The exponents of the eigenmode solution follow from the solution to det(C)=0 and 
show multiple branches as a function of 9, a, and r\. Asymptotically, the lowest branch 
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FIG. 8. Streamline plot of the eigenmode solutions (dark gray, blue online) and (light gray, 
red online) for decreasing a and increasing 9, for a viscosity ratio 77 = 1. 

that has i?e(A^) > 1 (to ensure regularity of the velocity field at the origin) dominates. 
The trivial roots A# = 1 and A# = 2 have to be excluded, as they do not represent an 
eigenmode, but are special cases for which the form (3) is degenerate. To understand the 
structure of the eigenvalues, it is insightful to first describe the limiting cases 77 = and 
77—7-00 when the eigenmode flows in liquids 1 and 2 are completely decoupled. In these 
two limiting cases, the leading-order exponents for the flows in both liquids are depicted as 
a function of 9 in Fig. 9(a) and (b) , for a given a = 90°. The leading-order exponent of 
the eigenmode flow in liquid 1 is shown as a black solid line and in liquid 2 as a gray solid 
line. In the decoupled case, the exponents for the eigenmode flows in liquids 1 and 2 are 
different. The expressions for the exponents are still relatively simple, see for example (21), 
compared to the case where the flows in liquids 1 and 2 are coupled. 

Introducing a weak coupling between the two liquids (dashed lines in Fig. 9(a) and (b)) 
ensures that the flows in both liquids obey the same exponent. This behavior means that 
the two separate branches represented by the solid lines in Fig. 9(a) and (b) give rise to two 
coupled eigenmodes. At the intersection point where the two separate branches (the solid 
lines) merge, a double root arises. Consequently, when the flow in the two liquids is coupled, 
there is a root-splitting behavior of the two coupled eigenmodes at this intersection point. 
A similar behavior was observed by Anderson & Davis 24 for a flow of two liquids between 
two rigid walls. The dominant eigenmode is the one that is closest to the lowest root of the 
two uncoupled solutions. This leading-order solution is shown as a dashed line in Fig. 9. 

To show how the viscosity ratio rj influences the exponent of the eigenmode flow in the 
coupled case, we plot the exponent of the eigenmode flow as a function of 9 for different 
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values of rj, for a given a = 90° in Fig. 9(c). The black solid line represents the leading-order 
exponent of the coupled eigenmode flow for rj = 1/100, which is the same case as depicted in 
Fig. 9(c) as the black dashed line. The opposite case, rj = 100, is shown as a black dashed 
line in both Fig. 9(b) and (c). We observe that as the viscosity ratio t] is increased for 
rj = 1/100 to r\ — 100, the maximum in the exponent becomes a minimum and vice versa. 

D. Dominant contribution 

The total evaporation-driven flow in the wedge geometry consists of three contributions: 
the flux, the moving contact line, and the eigenmode solutions. Each of the three contri- 
butions to the total flow in the wedge has a different scaling behavior with distance to the 
contact line: all are of the type p m , but with different m. In Fig. 10 these exponents are 
plotted as a function of angle a. The moving contact-line solution has exponent m = 1, 
which is the smallest of the three exponents for a < 90°. Therefore, the moving contact-line 
solution dominates the flow near the contact line for a < 90°. The exponent for the flux 
solution (in radians) is m = A = n/2a. Hence, for a > 90°, A < 1 and the flux solution 
is dominant. The leading-order exponent of the eigenmode solution varies not only with a, 
but also with 9 and r/. Figure 10 shows the result for 9 = 120° and 77 = 1. However, for all 
9, a, t] the eigenmode solution is always subdominant, since m — Xe > 1. 

IV. DISCUSSION 

We have derived self-similar solutions to the velocity field near the contact line of a liquid 
drop that evaporates from a liquid subphase. The nature of the flow strongly depends 
on the contact angles given by the Neumann construction. In both phases reversing flow 
structures can be observed, in particular when the angle a < 90°. We found that there 
are three contributions to the flow in the wedge geometry: one contribution that is driven 
by the evaporative flux from the drop surface, one contribution from the receding contact 
line, and one eigenmode contribution. The amplitudes of the flux and moving contact-line 
solution are known from the full spherical-cap solution for the vapor transport 13 . However, 
the amplitude of the eigenmode flow can only be obtained by solving the velocity field in 
the entire drop, which has to be done numerically. Which of the three flow contributions 
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FIG. 9. (a) Leading-order exponent of the eigenmode solution for liquid 1 (black, solid line) 
and 2 (gray, solid line) as a function of for r/ = 0, which means the flow in the two liquids is 
completely decoupled. The inset shows the boundary conditions experienced by liquid 1 and liquid 
2 in this case. The black dashed line represents the weakly coupled case where the flow in both 
liquids obeys the same exponent, for rj = 1/100. (b) Leading-order exponent of the eigenmode 
solutions for the decoupled case where r\ = oo. The black dashed line is the weakly coupled case 
for rj = 100. (c) Leading-order exponent of the coupled eigenmode solution for rj = 1/100 (black, 
solid line), rj = 1/10 (gray, solid line), 77 = 1 (light gray, dotted line), 77 = 10 (gray, dashed line), 
and rj = 100 (black, dashed line). For all graphs a = 90°. 



is asymptotically dominant close to the contact line depends on the contact angle a: for 
a < 90° we found that the moving contact- line solution dominates, whereas for a > 90°, 
the flux solution is dominant. The eigenmode solution is always subdominant, in contrast 
to the drop on a solid substrate (see 21 ), for which the eigenmode solution is the dominant 
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FIG. 10. Plot of the exponents of the flux (solid line), moving contact line (dashed line), and 
eigenmode (dotted line) solutions as a function of a. For a < 90° the moving contact-line solution, 
for which m = 1, dominates, for a > 90° the flux solution with m = A = ir/2a (in radians) is 
dominant. The leading-order exponent of the eigenmode solution, here plotted for 8 = 120° and 
ij = 1, is always subdominant (m = \e > 1). 

flow for contact angles larger than 133.4°. 

A fourth contribution to the flow could come from the Marangoni effect: the non-uniform 
evaporation of the drop leads to temperature gradients on the drop surface, which induces 
Marangoni stresses 5 . These Marangoni stresses can give rise to additional flow circulation 
in the drop, and also alter the velocity field in the liquid subphase. In principle, our analysis 
could be extended to include the Marangoni-driven contribution to the flow in the two liquids 
following the method of Ristenpart et al. 5 

It would be interesting to see if the reversing flow structures found in the present analysis 
can be measured in an experiment. A direct consequence of our analysis of the moving 
contact-line solution is the exact expression (including the prefactor) for the time evolution 
of the drop radius in time (20), which could easily be verified. To our knowledge, the 
evaporation time of a drop on a liquid subphase has not been investigated experimentally 
so far. 

We have applied the solutions to the flow in a wedge geometry that contains two liquid 
phases to the situation of a volatile drop evaporating from a non-volatile liquid subphase. 
The same solutions also hold for the inverse problem of a non-volatile drop on a volatile 
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subphase, such as oil drops on water. In that case, however, the outer problem for the 
evaporation is different, since we have a volatile liquid bath instead of a volatile drop. 
Evaporation will lead to a decrease in the bath level, leaving the drop geometry unchanged. 
Hence, the prefactor in the expression for the evaporative flux (8) will be different. Since 
there is no contact-line motion in this case, also the eigenmode flow can become important. 
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